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Graphane is obtained by perfectly hydrogenating graphene. There exists an intermediate material, partially 
hydrogenated graphene (which we call hydro graphene), interpolating from pure graphene to pure graphane. It 
has various intriguing electronic and magnetic properties . We investigate a metal-insulator transition, employing 
a quantum- site percolation model together with a graph theoretical approach. Hydrographene is an exceptional 
case in which electronic properties cannot be determined solely by the density of states at the Fermi energy. 
Though there are plenty of zero energy state in wide range of hydrogenation density, most of them are insulating 
states. We also demonstrate that it shows a bulk ferromagnetic property based on the Lieb theory. 



Graphene is one of the main topics in condensed matter 
physics because of its unusual electronic properties |Q] 0]. It 
is a promising material to design future nanoelectronic de- 
vices. To make a device it is necessary to confine electrons 
within a certain finite domain. However, this is impossible in a 
graphene sheet by applying electric field due to the Klein tun- 
neling effect, which is an intrinsic nature of Dirac electrons. 
One method is to cut graphene to form a required domain. Re- 
sultant graphene derivatives such as graphene nanoribbonsj3|] 
and graphene nanodisks[4] show remarkable electronic and 
magnetic properties. 

Recently graphane has been attracting much attention 
\$\. It is a graphene derivative obtained by perfectly hydro- 
genating graphene (FigO}. Graphene is a semimetal with 
each carbon forming an sp 2 orbital, while graphane is an in- 
sulator with each carbon forming an sp 3 orbital. Graphane 
provides us with an alternative method of confining electrons 
within a finite domain, since 7r-electrons are excluded from 
hydrogenated carbons. Namely, by dehydrogenating a finite 
region of a graphane sheet, electron can be confined within 
this region. By dehydrogenating a graphane sheet accord- 
ing to a circuit pattern, in principle it is possible to fabricate 
a graphene-graphane complex, embodying various electronic 
functions fl. 

In this paper, we explore an intermediate material, par- 
tially hydrogenated graphene, which interpolates from pure 
graphene to pure graphane H. Let us call it hydrographene. 
There are two types of hydrogenation, single-sided percola- 
tion and double- sided percolation. Employing a graph theo- 
retical approach to the site-percolation model 111 Oil , we present 
an intuitive and physical picture revealing the electronic and 
magnetic properties of hydrographene with hydrogenation 
density q, < q < 1. As q increases, there occurs a 
metal-insulator transition at q c . It is mapped to a ferromagnet- 
paramagnet transition in a Potts model with q c corresponding 
to a critical temperature T c . We remark that hydrographene is 
an exceptional case in which electronic properties cannot be 
determined solely by the density of states (DOS) at the Fermi 
energy. This is because there emerge a number of insulating 
states which contribute to the DOS at the Fermi energy but not 
to the conductivity. We also show that it is a bulk ferromagnet 
in the context of the Lieb theory. We find that single- sided hy- 
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FIG. 1. (Color online) (a) Illustration of graphene. There are carbons 
on a honeycomb lattice. They are grouped in two inequivalent sites 
A and B. (b) Illustration of graphane. Hydrogens are attached up- 
wardly to A sites, and downwardly to B sites. A honeycomb lattice 
is distorted. 



drogenation is more efficient to form large magnetic moment 
than double- sided hydrogenation. 

Model: In hydrographene a finite density of hydrogens are 
absorbed randomly to graphene[7]. Carbons with absorbed 
hydrogen form an sp 3 bond, where no 7r-electron exists. To 
simulate this fact we model the system by the Hamiltonian, 



(i) 



(id) 



where Cj (q) is the annihilation (creation) operator of n- 
electron, t = 2.7 [eV] is the transfer energy, and denotes 
the next-nearest neighbor sites of the honeycomb lattice. To 
exclude the electron at the hydrogenated site r we have im- 
posed an infinitely large on-site potential V at r. Thus the 
sum J2' r mns over the lattice sites where carbons are hydro- 
genated. 

We analyze a quantum site-percolation problem based on 
this Hamiltonian. Electromagnetic properties depend on per- 
colation networks (Figj2]), where hydrogenated carbons have 
been removed from the network. The resultant carbon net- 
work is nothing but a graph with vertices corresponding to car- 
bons. The model Hamiltonian has one-to-one correspondence 
to the adjacency matrix of a graph. Though the Hamiltonian is 
noninteracting, the problem is highly nontrivial because per- 
colation network is highly nontrivial. 
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FIG. 2. (Color online) Graph representation of percolation network 
in a honeycomb lattice, (a) 30% hydrogenation on the two sides. 
The system is not connected, (b) 30% hydrogenation on one side. 
The system is connected. 



Graphene is a bipartite system made of carbons at two in- 
equivalent sites (A and B sites) on a honeycomb lattice, as 
illustrated in Fig [Ha). These carbons can absorb hydrogens 
randomly. The reaction with hydrogen is reversible, so that 
the original metallic state and the lattice spacing can be re- 
stored by annealing [6D. However, the way of attachment is 
opposite between the A and B sites in graphane [FigQIb)]. As 
a matter of convenience, hydrogens are absorbed upwardly to 
A sites and downwardly to B sites. We can investigate two 
types of percolation problems in hydrographene. In one case, 
hydrogenation occurs randomly on both A and B sites. We 
call it double-sided percolation. In the other case, hydrogena- 
tion occurs only on A sites. We call it single- sided perco- 
lation. Single- sided hydrogenation can be manufactured by 
resting graphene on a Silica surface H. 

The number of lattice sites with no defects is N c . We as- 
sume N c is even for simplicity. We define the hydrogena- 
tion density by q — M/N c for double-sided hydrogenation 
and q = 2M/N C for single-sided hydrogenation, where M is 
the number of hydrogenated carbons, M = Y^ r - We choose 
the positions of hydrogenated carbons by the Monte Carlo 
method. Physical quantities are to be determined by taking 
the statistical average. 

Isolated carbons: We refer to carbons with no adjacent car- 
bons as isolated carbons, to those with one adjacent carbon 
as edge carbons, to those with two adjacent carbons as cor- 
ner carbons, and to those with three adjacent carbons as bulk 
carbons. We apply a combinatorial theory and estimate the 
numbers of these different types of carbons by using the fact 
that each site is occupied with probability q. We also use the 
quantity p = 1 — q in what follows. 

In the case of double-sided percolation, we calculate the 
number of isolated carbons to be Co = pq 3 N c , that of edge 
carbons to be C\ = 3p 2 q 2 N c , that of corner carbons to be 
C2 = 3p 3 qN c and that of bulk carbons to be C 3 = p 4 N c . 
They satisfy the relation Co + C\ + C 2 + C 3 = pN c . The total 
bond number is \ (d + 2C 2 + 3C 3 ) = %p 2 N c . 

In the case of single-sided percolation, we calculate the 
number of isolated carbons to be Co = ^q 3 N c , that of 
edge carbons to be C\ — ^pq 2 N c , that of corner car- 
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FIG. 3. (Color online) The numbers of isolated carbons C3, edge 
carbons C2, corner carbons C\ and bulk carbons Co in unit of N c . 
The horizontal axis is the hydrogenation density q. 
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FIG. 4. (Color online) The number of the zero energy states. The 
unit cell size is 100. The average is taken 1000 times. 



bons to be C2 = ^p 2 qN c and that of bulk carbons to 
be C 3 = \ (p 3 +p)N c . They satisfy the relation Co + 
C\ + C 2 + C 3 = \ (p + 1) iV c . The total bond number is 
\ (Ci + 2C 2 + 3C 3 ) = 3pN c . 

We show the numbers of various types of carbons estimated 
in this way in FigEl In order to verify that the combinatorial 
theory yields correct results, we have also carried out numer- 
ical calculations, and found that the results agree completely 
between them as in Figj3l 

Zero-energy states: We investigate the number of the zero- 
energy state A^o. It is given by diagonalizing the Hamiltonian, 
or equivalently by the difference between the dimension and 
the rank of the Hamiltonian, No = dim H—rankH. The anal- 
ysis is quite different between double- sided and single- sided 
percolations. 

We first discuss the double-sided case. We show the nu- 
merical results in Fig|4]a), where we have adopted the peri- 
odic boundary condition to the unit cell with sizes 100. We 
have also given the number of isolated carbons in the same 
figure. Note that one isolated carbon yields one zero mode. It 
is observed that most part of the zero modes have arisen from 
isolated carbons for high density hydrogenation limit (q « 1). 
A comment is in order. A metallic conductivity is propor- 
tional to Nq. However, this is not the case in hydrographene 
since isolated carbons cannot carry electric charges. We ex- 
pect a metal-insulater transition to occur, about which we give 
a detailed discussion soon after in this paper. 

We give physical interpretations of the zero modes. In the 
low density hydrogenation limit (q « 0), all sites are con- 
nected, i. e. the cluster number is only one. We can apply the 
Lieb theorem, which states that the number of the zero energy 
states is determined by the difference of the A- site and B-site. 
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FIG. 5. (Color online) The maximum cluster size (black dots) and the 
second maximum cluster size (red curves) for various hydrogenation 
density q in unit of N c . The unit cell size is 10000. The average is 
taken 100 times. 




(b) single-sided percolation 
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FIG. 6. (Color online) The mean finite cluster size for various hydro- 
genation density q. The unit cell size is 40000. The average is taken 
100 times. 



It is exactly obtained as 

M 



=0 



N c 



M — s ) \ M 



(2) 



in terms of binomial coefficients. The distribution is sym- 
metric at M = N c /2. We have shown Nq by the choice of 
N c = 100 in Figffia). By taking 7V C ->- oo, we obtain 



7V L = y/2N c pq/ir. 



(3) 



On the other hand, in the high density hydrogenation limit 
(q « 1), we can apply the dilute gas model. Every sites are 
isolated or unconnected with each other and give zero energy 
states. Then the number of zero is simply given by 



N c (l ~ q). 



(4) 



We have shown Nq by the choice of N c = 100 in Fig0|a). 
A simple interpolating formula reads as 



(5) 



where Co, Ci, C2 and C3 are the numbers of isolated carbons, 
edge carbons, corner carbons and bulk carbons, respectively. 
It explains the numerical result reasonably well, as in Figj4{a). 
The first and second terms represent the contributions from 
the isolated and edge carbons, respectively. The third term 
follows from the Lieb theorem with an appropriate correction. 

We now discuss the zero-energy states in single-sided per- 
colation. The number of the zero-energy states in each clus- 
ter is given by jVg luster = | JV|J uster - jv^ luster | in general. Here, 
jyduster > jyduster smce on i v a sites are hydrogenated. Hence, 
the total number of the zero-energy states is given by 



\ A jyduster 







N A = qN c . 



(6) 



This is confirmed perfectly by the numerical result in FigJJb). 
Thus, the numbers of the zero-energy states are very different 
between one-sided and double- sided hydrogenations. 

Metal-insulator transition: The bulk electronic properties 
of hydrographene is closely related to the connectivity of 
graph. We have calculated numerically the maximum cluster 



size N mc as a function of hydrogenation in the system with 
iV c = 10 4 , which we show in FigB The ratio N mc /N c is 
approximately equals to the percolation probability P(q). It 
is by definition the probability for an infinitely large cluster 
to appear in an infinitely large system. In the present system 
it may be interpreted as a probability that the right-hand side 
and the left-hand side of a hydrographene sheet is connected 
by a cluster. 

It is seen in Fig \5\ that the maximum cluster size decreases 
almost linearly as q increases, and suddenly becomes zero at 
a critical density q c . This is a characteristic feature of a phase 
transition. Indeed, it is a phase transition since the present 
percolation system is mapped to the ferromagnet as we point 
out later: See (|9]). For double-sided percolation the critical 
density is determined as q c = 0.32, which is consistent with 
the well-known result! 10] on percolation in the honeycomb 
lattice. For single- sided percolation the critical density is de- 
termined as q c = 0.55. 

The percolation probability is closely related to the conduc- 
tance of hydrographene. For q < q c , there is a large cluster 
around the sample. Thus the hydrographene is metallic. On 
the other hand, the hydrographene is insulator for q > q c be- 
cause there is no such cluster. Namely, a percolation-induced 
metal-insulator transition occurs at q = q c . It is intriguing that 
there is a large amount of the zero-energy states even though 
q > q c . This is highly contrasted with the normal metal, 
where the metal or insulator can be determined by the exis- 
tence of the zero-energy states. In hydrographene, the con- 
nectivity of carbon network plays a crucial role for the metal- 
insulator transition. 

We have also calculated numerically the mean size of finite 
cluster S (q) as a function of hydrogenation, which we show 
in FigO It diverge at q c = 0.31, indicating a critical behavior 



S(p)cx\q 



(7) 



for double- sided hydrogenation. 

Magnetization: We proceed to analyze the magnetic prop- 
erty of hydrographene. We are able to make a general argu- 
ment with respect to the magnetic moment. According to the 
Lieb theorem valid to the bipartite system at half-filling in 
connected graph, the magnetization M is determined by the 
difference between the numbers Na and Nb of the A and 
B sites in each cluster, M = \ \ Na — Nb | , where the spin 
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mapping The hydrogenation parameter q corresponds to 
the temperature T of a ferromagnet, 



FIG. 7. (Color online) The magnetization in unit of h/ 2 in the pres- 
ence (red dots) and absence (blue dots) of ferromagnetic coupling 
for various hydrogenation density q. A black curve represents the 
formula ([8}. Inset: An extended figure of the magnetization without 
ferromagnetic coupling. The unit cell size is 10000. The average is 
taken 500 times. 



direction is arbitrary in general. The magnetization of hydro- 
graphene is determined by that of the largest cluster, since it 
is dominant: See FigO We show the numerically calculated 
magnetization as a function of q in Fig (7] It is well approxi- 
mated by the relation 

M(q) = P(q)Nk, (8) 

where Nq is given by eq.® for double-sided hydrogenation 
and eq.© for single-sided hydrogenation: See Fig|7{b). The 
system is ferromagnet for q < q c , and paramagnet for q > q c . 

The maximum magnetization per site is obtained M — 
3 x 10 -3 at q = 0.21 for double- sided hydrogenation and 
M = 0.17 at q = 0.44 for single- sided hydrogenation. 
Single- sided hydrogenation is about 500 times more efficient 
than double- sided hydrogenation because only A sites are hy- 
drogenated in single- sided hydrogenation. The magnetization 
is proportional to the number of the sites. In other words, hy- 
drographene shows the bulk magnetism, which is highly con- 
trasted with the edge magnetism in graphene nanoribbons and 
nanodisks, where magnetization is proportional to the number 
of carbons along the edge. Bulk magnetism is desirable be- 
cause edge magnetism disappears in the thermodynamic limit. 

A comment is in order. It is intriguing that the magne- 
tization is zero for the perfect single- sided hydrographene 
(q = 1), which is in contradiction with the previous result flllll 
obtained based on a density functional theory study. This is 
not surprising since the direction of magnetization is arbi- 
trary in each cluster in our noninteracting model, as we have 
noted. The spin directions of two adjacent clusters may be 
aligned due to an exchange interaction if it is present, and it 
is probable that all clusters have the same spin direction. The 
maximum magnetization per site is obtained M = 0.13 at 
q = 0.68 for double-sided hydrogenation: See Fig|7]a). On 
the other hand, the magnetization increases linearly as a func- 
tion of the hydrogenation parameter q for single- sided hydro- 
genation: See FigHb). This reproduces the result[llj for the 
perfect single- sided hydrographene (q = 1). 

Percolation and ferromagnet: There exists a close rela- 
tion between the percolation problem and the ferromagnet. 
Indeed, the percolation problem is mapped to the zero- state 
Potts model on Kagome lattice via the Kasteleyn-Fortuin 



q = exp(-J/kT), 



(9) 



with J being the exchange stiffness. The percolation prob- 
ability P (q) is mapped to the magnetization as a function 
of temperature T. Thus, the metal-insulator transition cor- 
responds to the ferromagnet-paramagnet transition. Further- 
more, the critical behavior © corresponds to the Curie- Weiss 
law of the susceptibility. We anticipate various properties fa- 
miliar in ferromagnet to be translated into hydrographene via 
the Kasteleyn-Fortuin mapping. 

In summary, we have discovered various intriguing elec- 
tronic and magnetic properties of hydrographene based on 
a quantum-percolation model. We have proposed a new- 
type of percolation, single-sided percolation, which is qual- 
itatively different from the usual percolation on honeycomb 
lattice. Hydrographene is an ideal system to investigate from 
the viewpoint of quantum site-percolation transition. 
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